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Abstract 

We propose a spin cell based on photon-assisted tunneling through a conventional semiconductor 
barrier. The Dresselhaus spin-orbit interaction is included to break the spin rotation symmetry. 
Due to the in-plane electric field induced asymmetric momentum distribution in one lead, contin- 
uous flows of spin currents are driven through a barrier by a AC field. The net charge current 
remains zero. The spin current via photon-assisted tunneling can be readily adjusted via tuning 
the AC frequency or the in-plane electric field. This device may function as an ideal spin cell to 
supply spin currents in the spintronics circuit. 
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In recent years, an emerging technology named spintronics has attracted much attention 
(for review, see JlU3]). The motion of spin-polarized electrons produces a spin-dependent 
particle current y] . To form a complete spintronics circuit, it is vital important to have a 
spin cell which supplies continuous spin current through the circuit. Recently, such a spin cell 
was proposed based on the AC-driven double-dot system in a spatially nonuniform external 
mae netic field |3j. This spin cell can provide a spin current without a charge current by care- 
fully adjusting the device parameters. However, this spin cell met with some challenges such 
as the nonuniform magnetic field. An efficient non-magnetic spin cell is thus extremely de- 
sirable. The spin-orbit interaction due to inversion asymmetry|2j provides us the possibility 
to construct an electrical driving spin cell based on the conventional semiconductors. 

In this work, we propose an electrical driving single barrier spin cell, where the k 3 Dres- 
selhaus spin-orbit interaction in a confined region is employed to break the spin rotation 
symmetry. A continuous flow of spin current is generated when the isotropic momentum 
distribution of electrons is broken by an in-plane electric field, while the net charge cur- 
rent remains zero. By applying an AC field, the spin current can be greatly enhanced via 
photon-assisted tunneling and the spin injection can be readily adjusted via tuning the AC 
frequencies or the in-plane electric field. A tunable electrical driven spin cell can be achieved. 

The model device is illustrated in Fig. 1. The device consists of a piece of zinc blende 
semiconductor potential barrier which is confined by two ideal (free of charge and spin inter- 
action) conducting leads. The central region and the two leads are connected via tunneling. 
Let the growth direction along x || [0, 0, 1] direction. The coordinate axes x, y, and z are set 
to be parallel to the cubic crystallographic axes [0,0,1], [1,0,0], and [0,1,0]. The wave vector 
of electrons is k = (A^kn), where kii = (k y ,k z ) is the in plane wave vector and k x is the 
wave vector along the growth direction. We assume a perfect match at the interfaces of the 
device, where kii is a good quantum number. The Hamiltonian can thus be diagonalized in 
the basis ky . The central barrier we are interested is an open system where electrons can flow 
back and forth between the barrier and the two leads. The dynamics of this open system can 
be well described by an effective non-Hermitian Hamiltonian |4]: Tic = He — iTi, where He 
is the Hamiltonian of the central barrier and the non-Hermitian term — includes all the 
boundary conditions such as the barrier-lead coupling. By taking the assumption that the 
incident electron's kinetic energy is much smaller than the height of the barrier potential 5^ 
the central Hamiltonian He of the zinc blende semiconductor is then given in the nearest 
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neighbor tight-binding approximation |6( by 

h c = Yl K(k||)4 k||CT c nk||CT -t(4 k||(7 c n+lik||(T + F.a) (i) 
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where Ki(k||) contains all the on site energies at x = na (n = 1,2, ...,N), a is the dis- 
crete lattice constant, t = h 2 /2m*a 2 is the nearest neighbor hopping matrix element, 7 
is the Dresselhaus constant of the chosen material, and c^ k||CT (c„ k||0 -) creates (annihilates) 
a spin a electron at site n with a given parallel wave vector kn. In the wide band ap- 
proximation and symmetric coupling, the non-Hermitian part of Hamiltonian has the form 
S = Skyer f c ik||cr c ik||o- + f cjvjj^CATkno-, where r is the energy independent level width function. 

For a high barrier height, the electron transmission probability is greatly reduced by the 
nonresonant tunneling process. We apply an AC field to achieve a tunable spin cell via 
photon-assisted tunneling!?], Suppose the junction is illuminated by a monochromatic 
AC field of frequency uj. The on-site energies of the barrier are modified by the AC field. The 
total Hamiltonian is written as H = He + H^c, where Haq = — J2na A™ C0& {^t)^ n \L\\a c nk.\\a- 

In this work, we generalize the Floquet-Green formalismjq] to the spin-involved transport 
problems. For a given in-plane wave vector kn, solution of the time-dependent Schrodinger 
equation in the Floquet ansatz takes the form \ijj m {t)) = exp[—ie m t/h]\(p m (t)), where e m is 
the quasienergy. The Floquet states \(p m (t)) can be expanded in Fourier expansion form as 
Wm(t)) = J2 P \ fm,p) 6xp[— iput],p E Z. Both the Floquet states and the quasienergies are 
determined by the eigenvalue equation J r \(p m {'t)) = ^m\Vm(t)) with a non-Hermitian operator 
T — Ji — ihj^ 3|. Since the eigenvalues are generally complex, we have to solve the adjoint 
equation .F*|0 m (i)) = ej n |0 m (t)) to form a complete biorthogonal basis. The propagation of 
an electron from site n! with spin a' and energy e to the state at site n and spin a under 
the absorption (p > 0) or emission {p < 0) |p| photons is described by the retarded Green 
function 

rv / £ n - T (nv\¥m, P+P >)(<Pm,p>\n'a') , , 

m,p' t t m p (IW 

After obtaining the Green functions, spin-involved transport properties can be described 
by the 2x2 p-photon assisted transmission matrices T P LR and T^ L , where the matrix element 
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TaBaa' denotes the time-averaged p-photon assisted transmission probability for electrons 
from the (3 lead to the a lead with the initial and final spin state being a 1 and a respectively. 
These transmission matrices can be easily obtained by generalizing the result in Ref. 8J] 
to the present spin dependent transport problem. The diagonal and non-diagonal elements 
characterize, respectively, the spin-conserved and the spin-flip transport probabilities. Our 
model has a symmetric design and a standing wave AC potential. The two transmission 
probabilities are equal since the electron has no way to distinguish the left lead or the 
right lead. We therefore omit the subscript of the transmission probability as T p (e, ky) for 
simplicity. 

We apply an external electric field F = [F cos 9, F sin 9) in the y-z plane in the left lead to 
create an anisotropic lateral momentum distribution. We assume the central barrier is much 
higher than the potential caused by the electric field. In the weak electric field limit, the 
tunneling correction to the distribution function in the lead can be neglected. The difference 
between the left and the right distribution function in the relaxation time approximation is 
given byQ, Ell 

fL-f R *eT tr .^(v v F), (3) 

?iki 

where r tr is the transport relaxation time for electrons and vy = is the in plane electron 
velocity. At low temperature (T — > 0), the derivative — can approximately be replaced 
by a simple 5— function. The spin current can be found from the time evolution of the 
spin-resolved occupation number operator and the spin current density j a (in unit of A/m 2 ) 
is given by 

per 

where the in plane wave vector ky = (k\\ cos rj, k\\ sin 77) (77 is the polar angle) and the Fermi 
wave vector kp are determined by the carrier density in the leads. 

We numerically investigated the spin cell model with realistic parameters. We choose 
the material of central junction to be GaSb which has a larger Dresselhaus constant 7 
than other III-V materials 5]. The material parameters are given by7 = 187 eV-A 3 and 
m* = 0.041mo, where tuq is the electron mass. The length of the central junction in our 
device is L = 5a = lOOA.The AC potential parameter is fixed at A m = O.lt. The Fermi 
energy of leads and the barrier on site energy are fixed at and At, respectively. The carrier 
density at two leads is set to be 2 x 10 17 cm -3 and the Fermi wave vector is estimated to be 



4 



1.8 x 10 8 m -1 . A relaxation time of 1 ps is used in our numerical calculations. 

Fig. 2 depicts the surface spin current density as a function of the driving field frequency 
for a given external in-plane electric field. The in-plane electric field is F = 1.6 x 10 3 
V/m with the polar angle 9 = tt/3. The left lead electron distribution in the (k y ,k z ) space 
is therefore no longer isotropic. From the symmetry properties of our Hamiltonian, we 
have T°"(e, k||) = T~ a (e, — kn). Therefore, the spin current densities satisfy the relation: 
2a = —j-a- The two spin currents have the same magnitude but opposite signs. Therefore, 
the net charge current is always zero. When the frequency of the AC field matches the 
energy difference between the Fermi energy and the resonant level of the central junction, 
the spin current reaches its maximum value due to photon assisted tunneling. When the AC 
frequency is far off the resonant frequencies, the AC pumped spin currents decrease. Due 
to the spin-orbit interaction, the degeneracy of the lowest resonant level is lifted. One can 
find two photon-assisted tunneling resonant spin current peaks in Fig. 2. By scanning the 
frequency across the two resonant frequencies, the spin currents will reverse their direction. 
Therefore, we can readily control the direction and the magnitude of the spin current of the 
proposed spin current source by tuning the AC field frequency. If we decouple our device 
with the external circuit, spin accumulation occurs at the edges of the device due to the 
balance of spin accumulation and spin relaxation. This spin-accumulation may be exper- 
imentally detected by optical techniques based on Kerr/Faraday rotation with submicron 
spatial resolution [l^ . 

The output spin current can also be tuned by the in-plane electric field. In Eq. (J2J, the 
spin current density has a trivial proportion relation with the amplitude of the electric field. 
However, the polar angle of the in-plane electric field provides us another effective means to 
control the spin currents. In Fig. 3, we show the numerical calculated surface spin current 
density as a function of the polar angle 9 of the in-plane electric field. One can see that the 
spin currents are modulated by 9 in a sinusoidal way. This sinusoidal modulation is caused 
by the general feature and symmetry of the field induced nonequilibrium distribution and 
the spin-orbit interaction. We note the spin quantization direction is along z axis. From 
Eq. (@J), one can find that j a is modulated by cos9. For 9 = nir, n 6 Z, i.e. the electric 
field is vertical to the z direction (in y direction), the spin currents are suppressed to zero. 
However, when the in-plane electric field deviates from y axis, a non- vanishing spin current 
will be generated. Especially when the electric field is along the z axis (9 = ^y^-ti), this 
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spin current reaches its maximum. 

In conclusion, we have studied an AC driven spin cell by generalizing the Floquet-Green 
function method. The k 3 Dresselhaus spin-orbit interaction due to bulk inversion asymmetry 
contributes to the spin- dependent effective mass correction and leads to the splitting of the 
resonant level. Driving the lead electron distribution out of equilibrium by an in-plane 
electric field, a continuous flow of spin current can be generated by the AC field while the 
net charge current remains zero. The direction and the magnitude of the generated spin 
current can be readily tuned by changing the AC frequency or the in-plane electric field 
direction. Besides all the basic characteristics of spin cell presented in Ref . jsj , our spin cell 
functions without external magnetic field or any ferromagnets. It can be operated in a fully 
electrical way. By integrating our proposed spin cell with spin information processors, one 
can form a complete spintronics circuit based only on conventional semiconductors. 

This work was supported by Korea Research Foundation, Grant No.(KRF-2003-070- 
C00020). 
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FIGURE CAPTIONS: 



FIG. 1: Sketch of the proposed spin current source. The material growth direction is 
along x axis. The left lead electrons are driven by an in-plane (y-z plane) electric field and 
the central junction is illuminated by an AC field. 

FIG. 2: The surface spin current densities as a function of the driven frequency (in unit 
of the hopping matrix term t = h 2 /2m*a 2 ) for a given external in-plane electric field. 

FIG. 3: Numerical calculated surface spin current densities versus the polar angle 9 of 
the in-plane electric field at a fixed frequency Hcu = 2.261 
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FIG. 1: by B.H. Wu, and Kang-Hun Ahn 




FIG. 2: by B.H. Wu, and Kang-Hun Ahn 
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